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Abstract 

O 

■ Possibilities to search for 2(3 decay of a/ (3 unstable nuclei, many of which have higher energy 

release than "conventional" {(3 stable) 2(3 candidates, are discussed. First experimental half- 
life limits on 2(3 decay of radioactive nuclides from U and Th families (trace contaminants of 
the CaW0 4 , CdW0 4 and Gd 2 Si0 5 scintillators) were established by reanalyzing the data of 
low-background measurements in the Solotvina Underground Laboratory with these detectors 
(1734 h with CaW0 4 , 13316 h with CdW0 4 , and 13949 h with Gd 2 Si0 5 crystals). 
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1 Introduction 

Recent observations of neutrino oscillations with atmospheric [1], solar [2], reactor [3] and 
accelerator [4] neutrinos manifest the non-zero neutrino mass and prove an existence of new 
physical effects beyond the Standard Model (SM) of electroweak theory. The discovery of 
the neutrino mass gives an extraordinary motivation for experimental searches for neutrinoless 
(Ov) double beta {2(3) decay of atomic nuclei, an exotic process (forbidden in the SM), which 
violates the lepton number on two units [5, 6, 7]. While oscillation experiments are sensitive 
only to the neutrino mass difference, the measured Qv2(3 decay rate can give an absolute scale 
of the effective Majorana neutrino mass and, consequently, could test different neutrino mixing 
models. 

Up-to-date this process still remains unobserved, and only half-life limits for 0u mode were 
established in direct experiments (see reviews [5, 6, 7] and compilation of 2(3 decay results [8]). 
The highest bounds were obtained for the following nuclides: Ty 2 > 10 21 yr for 48 Ca [9], 150 Nd 
[10], 160 Gd [11], 186 W [12]; Ty 2 > 10 22 yr for 82 Se [13], 100 Mo [14]; Tfy 2 > 10 23 yr for 116 Cd [15], 
128 Te, 130 Te [16], 136 Xe [17]; and Tffc > 10 25 yr for 76 Ge [18]. These results have already brought 
the most stringent restrictions on the values of the Majorana neutrino mass (m u )< 0.2-2 eV, 
the right-handed admixtures in the weak interaction rj « 10~ 8 , A ~ 10~ 6 , the neutrino-Majoron 
coupling constant gu w 10 -4 , and the i?-parity violating parameter of minimal supersymmetric 
standard model ~ 10~ 4 [5, 6, 7, 8]. 

However, to make the discovery of the Oz/2/9 decay indeed realistic, the present level of the 
experimental sensitivity should be enhanced to (m„) ~ 0.01 eV. It is the great challenge, and 
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several large-scale projects of the next generation experiments on 0u2(3 decay search were pro- 
posed during few last years, like CAMEO, EXO, GEM, GENIUS, Majorana, MOON and others 
(see, f.e., reviews [5, 6, 7]). These projects intend to use up to tons of super-low background 
detectors made of enriched 2(3 isotopes, and plan to perform measurements with them during 
at least ~10 years. However, there is a chance that even for so ambitious projects the required 
advancement of the experimental sensitivity could be not reached in virtue of the restrictions 
of the current technologies. Therefore, investigation of other nonstandard approaches and dis- 
cussion of some new and even unusual ideas could be fruitful for the future of the 2f3 decay 
research. 

In the process of 2f3~ (or 2f3 + ) decay (see Fig. la) two electrons (or positrons) are emitted 
simultaneously, thus, an initial nucleus (A, Z) is transformed to (A, Z+2) or to (A, Z—2), which 
is, in principle, possible if mass of initial nucleus M(A, Z) is larger than mass of final nucleus 
M(A, Z ± 2). Nevertheless, quite often, to surprise, the statement can be found in a literature 
that 2/5 decay can occur in only case if: (a) M(A,Z) > M(A,Z ± 2); (b) ordinary (3 decay 
of the initial nucleus (A, Z) to the intermediate nucleus (A, Z ± 1) is forbidden energetically 
(that is M(A, Z)< M(A, Z ± 1)), or such a (3 decay (A, Z) -» (A, Z ± 1) is suppressed by a 
large difference in spin between parent (A, Z) and intermediate (A, Z ± 1) nuclei. However, it 
is clear that demand (b) is not mandatory - in fact, it just specifies conditions, which makes 
2/3 decay study convenient, because in the case of a (3 unstable parent nucleus, it would be 
extremely difficult to distinguish 2(3 decays from the intensive f3 background. To this effect, 
until now only "conventional" 2(3 decay candidate nuclei, satisfying both demands (a) and (b), 
were studied in direct experiments [5, 6, 7, 8]. 

However, investigation of the 2/3 decay of initially unstable {(3 or/and a) nuclei (see Fig. lb, 
lc), despite its complications, could be interesting too [8]. With this aim, let us consider formula 
for the Oz/2/3 decay probability, restricting to the neutrino mass mechanism only (right-handed 
contributions are neglected) [19]: (T^) -1 = F 0u ■ |NME| 2 • (m u ) 2 . Here, F 0u is the phase space 
factor (proportional to the fifth power of energy release in 2f3 decay, Qpp), and NME is the 
corresponding nuclear matrix element. Thus, the bound on the effective neutrino mass, (m u ), 
which could be derived from the experimental half-life limit, T^, can be expressed as follows: 
lim (m u ) ~ |NME| _1 • {Qpp • T-^} -1 / 2 . Last equation means that sensitivity of 2(3 experiment 

to the neutrino mass bound (for the equal NME and Ty£ limit) is proportional to the Q^J 2 , 
hence, the larger Qpp value, the more stringent (m u ) restriction could be derived. For any of 
69 "conventional" 2f3 ± decay candidate nuclei [8] this 2f3 energy release does not exceed «4.3 
MeV, 1 while for some 2f3 ± nuclides, not satisfying demand (b), the Qpp values are up to ten 
times larger [20]. The energy releases [20] for all possible 2/3 ± candidates (~2300) are shown in 
Fig. 2 in comparison with those available for "conventional" 2/3 nuclides. 2 For example, Qpp 
for 19 B (or 22 C) is equal ~43 MeV, therefore its 0z/2/5 decay rate would be ^4xl0 6 times faster 
than that for 76 Ge with Qpp 2 MeV (supposing equal NME-s). 

Obviously, because of very hard problem in accumulating a large amount of fast-decaying 
parent nuclei, and due to enormous difficulties in discrimination of 2f3 decay events from inten- 
sive (3 background, no schemes of 2{3 experiments involving "unconventional" candidates have 

In fact, 2 nuc lides from that list, 48 Ca and 96 Zr, are potentially (3 decaying, however, up-to-date their (3 
decays were not observed being strongly suppressed because of a large change in spin. These two nuclides have 
the highest Q m values: 4.272 MeV for 48 Ca and 3.350 MeV for 96 Zr. Only two other "conventional" candidates 
have Qpfi greater than 3 MeV: 100 Mo (3.034 MeV) and 150 Nd (3.368 MeV) [20]. 

2 Thus, "unconventional" schemes highly increase the list of potentially 2/3 decaying nuclides for studies. 
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Figure 1: Different configurations of (A, Z) — (A, Z + 1) — (A, Z + 2) nuclei: (a) "conventional" 
2/3~ triplet when mass of intermediate nuclide (A, Z + 1) is larger than that of initial (A, Z) 
and final (A, Z + 2) nuclei, and thus ordinary (3~ decay of (A, Z) is forbidden; (b) and (c) 
"unconventional" configurations when 2(3 decay is one of the branches of (A, Z) decay. The 
similar picture can be drawn also for 2(3 + decay. 



been proposed till now. Nevertheless, two more or less realistic methods could be suggested for 
these investigations. 

(1) Use of artificial unstable nuclides produced with accelerators or reactors 3 , which then will 
be used for 2/3 decay quest. Current possibilities to create radioactive ion beams on accelerators 
and accumulate them in storage rings are on the level of 10 19 nuclei per year [21]. As concerning 
the reactor-produced isotopes, such an approach was already employed once in radiochemical 
search for 2(3 decay of 244 Pu where the limit (on all modes) was set as Ti/ 2 > 1.1 xlO 18 yr [22]. 

(2) Use of appropriate unstable nuclides in natural radioactive U/Th families. These chains 
are present in some quantities, as contamination, in any materials including those used for 
low background detectors. Thus, limits on 2(3 decays of such "unconventional" isotopes can be 
derived as by-product of any low background measurements with the proper detector (including 
searches for "conventional" 2(3 decay). Moreover, detectors specially loaded by these radioactive 
isotopes could be produced for such experiments in order to increase the number of nuclei for 
investigation. 

Further, we would like to concentrate on the second possibility and present the experimental 
results, which were obtained by using this method for the first time. 

3 Besides, such nuclides (but in much less quantities) could be produced in detector by cosmic rays. Nucleon- 
and muon-induced reactions in a target with given (A, Z) could result in creation of many different radioactive 
nuclei with lower A and Z numbers. 
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Figure 2: Energy release in double beta decay. The Qpp values were extracted from [20]. For 
"conventional" 2(3 decaying nuclides, points are shown in bold and connected by line. 



In order to make a choice of a right candidate for study, let us consider an experimental 
sensitivity of the mentioned approach. For radioactive chain in equilibrium, decay rates of 
different isotopes, R a/f3 = dN/dt, are the same. One can determine the number of nuclei N 
for each isotope by using the relation N = R a ^ ■ T™j 2 / In 2. Here a small contribution from 2(3 

decay is neglected; T™! 2 is the isotope's half-life for the usual a or (3 decay. From the other 

side, the limit on half-life for 2(3 decay is equal limT-j 2 ^ = e • In 2 • N • tj limS*, where e is the 
efficiency to detect the 2(3 process, t is the time of measurements, and lim S is the limit on the 
number of observed 2(3 events which can be excluded with a given confidence level on the basis 
of experimental data. Combining these two equations and expressing the a/ (3 decay rate R a ^ 
through the observed specific activity A a IP = R a ^ jm (m is the mass of detector), one can get 
finally: 

lim Ty 2 = e-m-t- A a/f3 ■ T?/£ / lim S. ( 1 ) 

From the latter one can see immediately that the larger is T^j^ half-life of candidate nuclide, 
the higher T^ 2 limit could be established. To obtain non- trivial bound lim 7^ > T"//, the 
condition e ■ m ■ t ■ A a ^ / \imS > 1 should be fulfilled. For typical values of e ~ 1 (if 2(3 
decaying isotope is containing in the detector itself), one year of measurements, limS' = 2.4 
counts (i.e., for zero observed 2(3 events [23]), m — 1 kg, and A a ^ = 1 /iBq/kg (so, for quite 
low contamination of detector by U/Th chains), we obtain the ratio \imTy 2 /Ty 2 ^ m 10. For 
the higher level of contamination A a ^ = 1 mBq/kg, this ratio is equal 10 4 . 
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In conclusion, to reach the higher sensitivity one has to select the candidate nuclides with 
the largest Qpp energies and with the longest T^j 2 values. In Table 1, where the properties 
of all potentially 2/3 decaying nuclei in U/Th families are summarized, there are 11 double (3 
unstable nuclides in the 232 Th chain, 6 nuclides in the 235 U chain, and 16 ones in the 238 U chain 
[24]. Unfortunately, most of these candidates are out of our interest because of their small Qpp 
energies, short T^j 2 half-lives, or because of tiny branching ratios in the chain, A. For example, 

210 T1 has one of the highest Q@p value 5.548 MeV, but its half-life is too short (Ty 2 = 1.30 m) 
and branching ratio is very low (A = 0.021%). 

In the present work we report the first experimental results of the quest for 2(3 decays 
of unstable nuclides in U/Th chains performed in the real-time measurements with the low- 
background CaW0 4 , CdW0 4 and Gd 2 Si0 5 crystal scintillators, which contain these nuclides 
as trace contaminations. 



2 Experiments and data analysis 

The experiments were carried out in the Solotvina Underground Laboratory of INR in a salt 
mine 430 m underground (~1000 mwe, with a cosmic muon flux of 1.7xl0~ 6 cm~ 2 s _1 , a 
neutron flux <2.7xl0~ 6 cm~ 2 s _1 , and a radon concentration in the air <30 Bq m -3 ) [26]. 
The total exposure was equal 1734 h with CaW0 4 (13316 h with CdW0 4 and 13949 h with 
Gd 2 SiOs) detector. 

Since detailed descriptions of the apparatus and technique of the experiments with CdW0 4 
and Gd 2 Si0 5 detectors are given in refs. [12, 15, 11], below we will describe only the main 
features, performances and experimental procedure used with the CaW0 4 crystal scintillator 
[27]. 

The CaW0 4 crystal used has dimensions 40x34x23 mm (mass of 189 g). Its measured 
light output (peak emission at 440 nm with decay time of ~9 /is) is «19% as compared with 
that of Nal(Tl). The crystal was viewed by the special low-radioactive 5" photomultiplier tube 
(EMI D724KFLB) through the quartz light-guide 010x33 cm. The detector was surrounded 
by a passive shield made of teflon (thickness of 3-5 cm), plexiglass (6-13 cm), high purity 
copper (3-6 cm), lead (15 cm) and polyethylene (8 cm). Two plastic scintillators (120x130x3 
cm) installed above the passive shield were used as a cosmic muons veto. Event-by-event data 
acquisition (the same used with CaW0 4 , CdW0 4 and Gd 2 Si0 5 detectors) records information 
on the amplitude (energy), arrival time and pulse shape of a signal (for Gd 2 SiOs pulse shape 
data were not recorded). For the latter, a transient digitizer based on the fast 12 bit ADC 
was used with the sample frequency of 20 MS/s [28]. Events in the chosen energy interval (for 
the CaW0 4 usually with the energies higher than «180 keV) were recorded in 2048 channels 
with 50 ns channel's width. The linearity of the energy scale and resolution of the detector 
were measured with 60 Co, 137 Cs, 207 Bi, 232 Th and 241 Am 7 sources in the energy range of 60- 
2615 keV. For example, FWHM is equal 9.7% at the energy 662 keV. The a/(3 ratio 4 of the 
crystal was measured with the help of collimated a particles ( 241 Am) in the range of 0.5-5.3 
MeV by using the set of thin mylar films (see for details [29]). Besides, a peaks of 147 Sm 
and nuclides from the 232 Th, 235 U and 238 U families, present as trace in the CaW0 4 crystal, 
were used to extend the energy interval up to ~8 MeV. Alpha peaks of 147 Sm, 210 Po, 218 Po, 

4 The a/ (3 ratio is defined as ratio of a peak position in the energy scale measured with 7 sources to the 
energy of a particles. 
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Table 1: List of candidate nuclides (present in U/Th radioactive chains), which could undergo 
2[3~ decay, and half-life limits, Ty 2 , determined in this work. The branching ratios, A, for each 
nuclide in chain are taken from [24]. Only for 238 U the half-life for 2/3 decay (all modes) was 
previously measured in radiochemical experiment as Ty 2 = (2.0±0.6)xl0 21 yr 



Parent 


Main channel of 


A, % 




limT 2 ^ at 68% C.L. 


nuclide 


decay and T™j 2 




MeV 


for Ou (2u) mode 


232 Th chain 


232 T > 
90 111 


a 


1.405 xlO 10 yr 


100 


0.842 


1.6xlO n (2.1xl0 9 ) yr^ 


8g Ra 


0- 


5.75 yr 


100 


2.173 




228 A 
89 AC 


P~+a 


6.15 h 


100 


0.016 




224p 


P~ 


3.30 m 


5.5xl0~ 6 


1.417 




220 Rn 
86 1X11 


a 


55.6 s 


100 


0.344 




2 ^Po 


a 


0.145 s 


100 


1.534 




2 iPb 


P~ 


10.64 h 


100 


2.828 


6.7 (0.4) yr**) 


2 iBi 


P~+a 


60.55 m 


100 


0.500 




™Pb 


j3~+a 


22.3 yr 


4.3xl0~ 9 


1.226 




208T1 
81 11 


P~ 


3.053 m 


35.94 


2.121 




206tt„ 
80 n & 


P~ 


8.15 m 


8.2xl0- 17 


2.841 




235 U chain 


90 - Ln 


P~ 


25.52 h 


100 


0.030 




2 8 23 Fr 


p~+a 


21.8 m 


1.380 


0.563 




2 ^5 9 At 


a+p~ 


56 s 


8.3xl0" 5 


1.918 




2 83 Bi 


P~ 


7.6 m 


8.0xl0~ 5 


2.971 




2 MPo 




1.781 ms 


100 


0.639 




2 ^Pb 


P~ 


36.1 m 


99.99977 


1.952 


1.5 (0.07) d 


238 U chain 


238 tt 
92 u 


a 


4.468 xlO 9 yr 


100 


1.147 


l.OxlO 12 (8.1xl0 10 ) yr 


234 Th 

90 111 


P~ 


24.10 d 


100 


2.470 


144 (2.5) yr 


234p„ 


P~ 


6.70 h 


100 


0.387 


17 (0.9) yr 

4.1xl0 4 (4.5xl0 3 ) yr 


2 |gRa 


a+ u 6 C 


1600 yr 


100 


0.476 


2 fjRn 


a 


3.8235 d 


100 


2.057 


2.8 (0.11) yr 


2 ^Po 


a+p- 


3.10 m 


100 


3.147 


3.7 (0.5) d 


2 ^ 5 8 At 


a+p- 


1.6 s 


0.020 


1.041 




2 tJpb 


P~ 


26.8 m 


99.980 


4.295 


87 (1.5) d 


2 ^Bi 


p~+a 


19.9 m 


100 


2.182 


29 (3.5) d 


2 iPb 


P~ 


10.64 h 


3.2xl0- 9 


2.828 




2 iBi 


p~+a 


60.55 m 


3.2xl0~ 9 


0.500 




210T1 

81 11 


P~ 


1.30 m 


0.021 


5.548 




2 ^pb 


P~+a 


22.3 yr 


100 


1.226 


1.2xl0 5 (8.6xl0 3 ) yr 


208T1 
81 11 


P~ 


3.053 m 


1.2X10" 9 


2.121 




2 8 °o 6 Hg 


P~ 


8.15 m 


5.6xl0- n 


2.841 






P~ 


3.38 m 


5.6xl0^ n 


7.986 




*) The limit was set with the 116 Cc 


1W0 4 ( ** } - 


with the Gd 2 Si0 5 ) detector. 
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232 Th, 238 U were selected with the help of a pulse-shape analysis (see below), while 214 Po, 215 Po, 
216 Po, 219 Rn and 220 Rn peaks were reconstructed with the help of a time-amplitude analysis 
from the data accumulated in the low-background measurements. All results (with external 
and internal a sources) are in a good agreement. In the energy interval 0.5-2 MeV the a//3 
ratio decreases with energy: a/ (3=0.21(3)— 0.020(15) E a , while it increases in the interval of 2-5 
MeV: a//3=0.129(12)+0.021(3)£ a , where E a is in MeV. 

Pulse-shape discrimination. The pulse shape of CaW0 4 scintillation signals can be 
described by formula: f(t) = J2 Mjiji— To) x (e~'/ ri — e~*/ r °), t > 0, where Ai are amplitudes 
(in %) and r, are decay constants for different light emission components, r is integration 
constant of electronics (~ 0.2 /j,s). The following values were obtained by fitting the average 
of a large number of individual pulses: Af=84%, t" = 9.1 /is, ^=16%, = 2.5 /is for «4.6 
MeV a particles and AJ=92%, t?=9.2 /is, A%=8%, r 2 7 =2.6 /is for wl MeV 7 quanta. This 
difference allows one to discriminate j(/3) events from those of a particles. For this purpose 
we used the method of the optimal digital filter (for the first time proposed in [30]), which 
previously was successfully applied with CdWC>4 scintillators [28]. To obtain the numerical 
characteristic of CaW0 4 signal, so called shape indicator (SI), each experimental pulse f(t) 
was processed with the following digital filter: SI = J2f(tk) x P(tk)/ f(tk), where the 
sum is over time channels k, starting from the origin of pulse and up to 75 /is, f(tk) is the 
digitized amplitude (at the time tk) of a given signal. The weight function P(t) is defined as: 
P(t) = {f a (t) - fj(t)}/{f a (t) + fj(t)}, where f a (t) and f 7 (t) are the reference pulse shapes 
for a particles and 7 quanta, resulting from the average of a large number of experimental pulse 
shapes. 

Shape indicator SI measured by the CaW0 4 crystal (40x34x23 mm) for alpha particles 
in the 1—5.3 MeV region does not depend on the direction of a irradiation relative to the 
crystal axes. Similarly, no dependence of the SI on 7 quanta energy (within the energy range 
0.1—2.6 MeV) was observed. The distributions of the shape indicator measured with a particles 
(E a :=» 4.6 MeV) and 7 quanta (~1 MeV) are depicted in the inset of Fig. 3 (the larger value 
of the shape indicator corresponds to the shorter decay time of scintillation pulse). As it is 
seen, distinct discrimination between a particles and 7 rays ((5 particles) was achieved. An 
illustration of the PS analysis of the background data, accumulated during 171 h with CaWC>4 
detector, is shown in Fig. 3 as scatter plot for the SI values versus energy. In this plot one can 
see two clearly separated populations: the a events, which belong to U/Th families, and j(/3) 
events. 

Background and radioactive contamination of the CaWC>4 crystal. As it is men- 
tioned earlier, the CaWC>4 crystal (40x34x23 mm) was measured with the help of the low 
background set-up, installed in the Solotvina Underground Laboratory. The energy resolution 
of the detector was determined with several 7 sources ( 60 Co, 137 Cs, 207 Bi, 232 Th and 241 Am) 
and can be fitted by function: FWHM 7 (keV) = -3 + ^6 .9Ej, where E y is energy of 7 quanta 
in keV. The routine calibrations were carried out weekly with 207 Bi and 232 Th sources. 

The energy spectrum of the CaWC>4 detector, measured during 1734 h in the low background 
apparatus was separated into a and (3 spectra with the help of the pulse-shape discrimination 
technique. 

First, the a spectrum, which is depicted in Fig. 4, has been analyzed. The intensive and 
clear peak at the energy 1.28 MeV (in 7 scale) can be attributed to intrinsic 210 Po (daughter 
of 210 Pb from 238 U family) with activity of 0.291(5) Bq/kg. Apparently, the equilibrium of 
the uranium family in the crystal was broken during crystal production, because peak of 238 U 
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Figure 3: Scatter plot of the shape indicator SI versus energy for 171 h background data 
measured with the CaWC>4 crystal scintillator (40x34x23 mm). Lines show ±2cr region of SI 
for 7 (/3) events. (Inset) The SI distributions measured in calibration runs with a particles 
(E a = 4.6 MeV which corresponds to ~1 MeV in 7 scale) and 7 quanta (~1 MeV). 



(see Inset (a) in Fig. 4) corresponds to the much lower activity of 14.0(5) mBq/kg. Peaks of 
the uranium's daughters 234 U, 230 Th, 226 Ra are not resolved (their Q a values are very close), 
however, the area of the total peak (at wl.l MeV) is in satisfactory agreement with the activity 
of 238 U and 226 Ra (see below the results of time-amplitude analysis for 226 Ra). Another member 
of the family, 222 Rn, is not discriminated from 210 Po (an expected energy of a peak is ~1.34 
MeV in 7 scale), while 218 Po is well resolved. 

In the low energy part of alpha spectrum (Inset (b) in Fig. 4) the peak at the energy «0.8 
MeV can be attributed to 232 Th with activity 0.69(10) mBq/kg. Weak alpha peak with the 
energy in 7 scale of 395(2) keV (corresponds to energy of a particles 2243(9) keV) can be 
explained by trace of 147 Sm (E a = 2247 keV, Ti/ 2 =1.06xlO n yr, isotopic abundance is 15.0% 
[31]) with activity 0.49(4) mBq/kg. The total internal alpha activity in the calcium tungstate 
crystal is ^0.4 Bq/kg. 

Besides, the raw background data were analyzed by the time-amplitude method, when the 
energy and arrival time of each event were used for selection of some decay chains in 232 Th, 
235 U and 238 TJ families. For instance, the following sequence of a decays from the 232 Th family 
was searched for and observed: 220 Rn (Q a = 6.40 MeV, T 1/2 = 55.6 s) -> 216 Po (Q a = 6.91 
MeV, T1/2 = 0.145 s) — > 212 Pb (which are in equilibrium with 228 Th). Because the energy of a 
particles from 220 Rn decay corresponds to ~1.6 MeV in 7 scale of CaWC>4 detector, the events 
in the energy region 1.4 - 2.2 MeV were used as triggers. Then all events (within 1.4 - 2.2 
MeV) following the triggers in the time interval 20 - 600 ms (containing 85.2% of 216 Po decays) 
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Figure 4: Energy spectrum of a events selected by the pulse-shape analysis from background 
data measured with the CaWC>4 detector (0.189 kg, 1734 h). (Inset a) The same spectrum but 
scaled up. It is well reproduced by the model, which includes a decays of nuclides from 232 Th 
and 238 U families. (Inset b) Low energy part of the a spectrum. 



were selected. The obtained a peaks (the a nature of events was confirmed by the pulse-shape 
analysis described above), as well as the distributions of the time intervals between events are 
in a good agreement with those expected for a particles of 220 Rn — > 216 Po — > 212 Pb chain 5 (see 
Fig. 5). On this basis the activity of 228 Th in the CaWC>4 crystal was calculated as 0.6(2) 
mBq/kg, which is in a good agreement with activity of 232 Th determined from a spectrum - 
0.69(10) mBq/kg. 

For the analysis of the 226 Ra chain ( 238 U family) the following sequence of (3 and a decays 
was used: 214 Bi (Q p = 3.27 MeV) -> 214 Po (Q a = 7.83 MeV, T 1/2 = 164 /is) -> 210 Pb. For 
the first event the lower energy threshold was set at 0.18 MeV, while for the second decay the 
energy window 1.6 — 2.4 MeV was chosen. Time interval of 90 — 500 /is (56.3% of 214 Po decays) 
was used. The obtained spectra for 214 Bi and 214 Po lead to the 226 Ra activity in the CaW0 4 
crystal equal to 5.6(5) mBq/kg. 

Finally, let us analyze the energy spectrum of /3( , y) events selected with the help of the 
pulse-shape technique and presented in Fig. 6. The counting rate for the /3(7) spectrum above 
the energy threshold of 0.2 MeV is f^O.45 counts/ (s-kg). The contribution of the external 7 rays 
to this background rate was estimated as small as ~2%, by using results of measurements with 
CdWC>4 crystal (mass of 0.448 kg) installed in the same low background set-up. Therefore, the 
remaining ^(7) events are caused by the intrinsic contaminations of the CaW0 4 crystal. Major 

5 The a peak with the energy E a «7.3 MeV, which is present in the energy distribution of the second event, 
can be attributed to 215 Po from 235 U family. Corresponding activity of 227 Ac in the crystal is 1.6(3) mBq/kg. 
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Figure 5: Alpha peaks of 220 Rn and 216 Po selected by the time-amplitude analysis from the data 
accumulated with the CaW0 4 detector. (Inset) The distribution of the time intervals between 
the first and second events (dots) together with the exponential fit (line). Obtained half-life of 
216 Po (0.17 ± 0.04 s) is in an agreement with the table value: 0.145(2) s [24]. 



part of this (3 activity can be ascribed to: 210 Bi (daughter of 210 Pb) - 0.291 Bq/kg; 234m Pa 
- 0.014 Bq/kg; 214 Pb and 214 Bi ( 238 U family) - ^0.007 Bq/kg; 211 Pb ( 235 U family) - ^0.002 
Bq/kg. The remaining (3 events can be explained by other (3 active impurities ( 40 K, 90 Sr in 
equilibrium with 90 Y, 137 Cs, etc.) probably present in the crystal. 

The summary of the measured radioactive contamination of the CaW0 4 , CdW0 4 and 
Gd 2 Si0 5 crystal scintillators (or limits on their activities) is given in Table 2 [11, 15, 27]. One 
can see that radioactive impurities in the CaW04 and Gd2Si05 crystals (available at present) 
are much higher (by factor of 10 - 10 3 ) than those of the CdW0 4 scintillators. 

3 Limits on 2/3 decay and conclusions 

The background spectra of (3/^ events selected by the pulse-shape discrimination method with 
the CaW0 4 , CdW0 4 and Gd 2 Si0 5 crystal scintillators were analyzed to search for 2/3 decays of 
unstable isotopes in U/Th chains. In general, we did not find any peculiarities in the measured 
spectra which can be attributed to the double beta decay processes searched for. Therefore, 
only the lower half-life limits were established in accordance with the equation (1). 

Values of lim S were simply estimated as a square root of the number of counts in the corre- 
sponding energy windows of the background spectra. For instance, in the spectrum measured 
with the CaW0 4 detector, there are 3937 counts in the energy interval 1110-1200 keV where 
the peak of the 0z/2/3 decay of 238 U is expected. Number of 238 U nuclei is N = 5.4xl0 14 , and 
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Figure 6: Energy spectrum of j3(j) events selected by the pulse-shape analysis from background 
data measured with the CaWCU detector (0.189 kg, 1734 h). It is described by (3 spectra of 
210 Bi (Qp = 1.16 MeV), 214 Pb (Qp = 1.02 MeV), 211 Pb (Qp = 1.37 MeV), 234m Pa (Qp = 2.27 
MeV), and 214 Bi (Qp = 3.27 MeV). (Inset) The energy region of 0z/2/3 decay of 238 U is shown 
together with expected peak of 0u2(3 decay which corresponds to Ty 2 = 10 11 yr. 



efficiency for this energy region is e = 0.85. The limS 1 = 63 counts together with the time of 
measurements t = 1734 h lead to the restriction on the half-life: T®% > 1.0 xlO 12 yr at 68% 
CL. 6 Considering the energy interval 200-760 keV (number of background counts is 4.7xl0 5 , 
efficiency e = 0.75), the half-life limit relatively to the two-neutrino 2(3 decay of 238 U was esti- 
mated as: Ty 2 > 8.1 xlO 10 yr at 68% CL. In the inset of Fig. 6, the energy spectrum of the 
CaW04 detector in the region of Qv2(3 decay of 238 U is shown together with the expected peak 
of 0i/2/3 decay (T 1/2 = 10 11 yr). 

Similarly, such a procedure was applied to estimate the T ly / 2 limits for 2(3 decays of other 
nuclides in U/Th radioactive families in the CaW04 crystal (from Table 1), as well as by 
using the data accumulated during 13316 h with the 116 CdW04 detector [12, 15] and with the 
Gd2Si05 scintillator (13949 h) [11]. Pulse-shape and time-amplitude analysis of events were also 
used to reduce backgrounds. For example, in 2(3 decay of 214 Pb (Q 2 ,g=4.3 MeV) the daughter 
nucleus 214 Po quickly (Ti/ 2 =164 /is) a decays with Q a =7.8 MeV. Such a chain of two events, 
first of which (with the energy 4.3 MeV) has the shape of scintillation flash specific for j/{3 
particles, and the second one (7.8 MeV, occurred within «1 ms after the first) has the pulse 



6 Similar results were obtained in more sophisticated approach, when the experimental spectrum was fitted 
by the sum of background model and response function of the set-ups for 2/3 decay searched for (both were 
Monte Carlo simulated with the GEANT3.21 and DECAY4 codes [11, 12, 15, 29]). Therefore, notwithstanding 
its simplicity, the "square root" approach gives a right scale of the sensitivity of an experiment. 
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Table 2: Radioactive contaminations measured in the CaWO"4 (mass of 189 g, measuring time 
1734 h), CdW0 4 (330 g, 13316 h) and Gd 2 Si0 5 (635 g, 13949 h) crystal scintillators. 



Chain 


Sub-chain 


Activity (mBq/kg) 








CaWC-4 


lie CdW0 4 


Gd 2 Si0 5 


232 Th 


232 Th 

228 Ra . 


.. 228 Ac 


0.69(10) 


0.053(9) 
< 0.004 


< 6.5 
< 9 




228 Th . 


.. end of chain 


0.6(2) 


0.039(2) 


2.287(13) 


235 JJ 


23l Pa 








< 0.08 




227 Ac . 


.. end of chain 


1.6(3) 


0.0014(9) 


0.948(9) 


238 JJ 


238 TJ 


234 Pa 


14.0(5) 


< 0.6 


< 2 




234 u 






< 0.6 






230 Th 






< 0.5 


< 9 




226 Ra . 


.. 214 Po 


5.6(5) 


< 0.004 


0.271(4) 




21 °pb . 


. . end of chain 


291(5) 


< 0.4 


< 0.8 



shape character to a's, gives a very distinctive signature, which helps to suppress background 
greatly. 

All results of the search for the 2(3 decays of the radioactive nuclides in thorium and uranium 
families obtained in this work are summarized in Table 1 (only half-life limits are given whose 
values exceed one day). 

In conclusion, we have presented results of the first experimental search (in direct counting 
experiment) for the 2/3 decays of a/ (3 decaying nuclides in U/Th chains, which are present in 
the CaW04, CdW04 and Gd2SiOs crystal scintillators as contaminations. The results were ob- 
tained by reanalyzing the data accumulated during few years of measurements in the Solotvina 
Underground Laboratory with these detectors, whose properties (radioactive contaminations, 
energy resolutions, a/ (3 ratios and pulse-shapes) were carefully studied. 

The half-life limits derived in this first attempt are in the range from a few days to 10 12 
yr, i.e. much more modest than those of "conventional" 2(3 decay experiments. Obviously, the 
sensitivity of our method could be further enhanced by producing scintillators loaded by thorium 
or uranium (the level of allowed concentration of these nuclides is restricted by the requirement 
of reasonable counting rate and by demand to keep satisfactory scintillation characteristics of 
the detector). For example, supposing the use of some fast scintillator (with decay time in 
the range of ns), the activity of the U/Th admixture could be increased to ~10 4 Bq/kg in 
comparison with the current level of ~1 mBq/kg. Together with the total mass of detectors 
enlarged from ~0.1 kg to 100 kg, it could allow one to advance the current limits on 2(3 half-lives 
of nuclides in U/Th chains by ~10 orders of magnitude, which looks interesting. 

Let us consider also the unstable nuclides, which can be produced with accelerators or 
reactor for the 2(3 decay searches. There exist plans of large-scale experiments with quickly (3 
decaying radioactive ions [21] (producing intensive beams of pure v e and u e ) for high precision 
measurements of neutrino oscillations with the future megaton Frejus detector [34]. These 
experiments could give, as by-product, some results on Oz/2/3 decay of involved nuclei. However, 
even principal schemes for extraction of such by-products were not debated yet. As an inspiring 
example, we can mention the experiment [35], where a active 221 Ra isotope with Ti/ 2 = 28 s 



12 



was produced by bombarding a thorium target with the 600 MeV proton beam. Then, the 
cluster radioactivity (emission of 14 C) of 221 Ra was observed with 7\/ 2 = 7.8xl0 5 yr [35], that 
is ^10 12 times longer than half-life of 221 Ra alpha decay. 

Besides, by analyzing the table of isotopes [24] for the long lived unstable nuclei (with T ± / 2 
> 1 yr), which can simultaneously undergo 2(3 decay, three nuclides were found with Q 2 p values 
higher than 4 MeV: 42 Ar (T 1/2 =32.9 yr, Q 2/3 =4125 keV), 126 Sn (T 1/2 ~lxl0 5 yr, Q 2/3 =4050 keV) 
and 208 Po (71/2=2.9 yr, Q 2EC =4280 keV). From them, 126 Sn seems to be the most interesting 
due to its longest life-time. As a possible detector for the 0u2(3 decay quest of 126 Sn, the liquid 
scintillator loaded by 126 Sn (up to ~10-20% in mass [36]) could be considered. The contribution 
of two successive single (3 decays ^Sn (T 1/2 ~lxl0 5 yr, Q =38O keV) -> 12 ^Sb (T 1/2 =12.5 d, 
Q/3=3670 keV) — > 12 |Te to background in the energy window of the 0z/2/3 decay of 126 Sn would 
be negligible due to the long half-life of the intermediate 126 Sb nucleus and taking into account 
the energy distributions of both ( 126 Sn and 126 Sb) (3 spectra. We recall that there is only one 
"convenient" 2(3 candidate with Q 2 p > 4 MeV: 48 Ca, whose natural abundance is very low 
(0.19%) that results in the world reserve of ~50 g of enriched 48 Ca isotope available for the 
investigations [9]. Note that for equal experimental 7\/ 2 limits on 0z/2/3 decay, an experiment 
with 126 Sn would be more sensitive to the neutrino mass than that with 48 Ca due to higher 
phase space available (which depends not only on the Q 2 p but on the Z value of a nucleus as 
well). However, perspectives to obtain considerable amounts of 126 Sn isotope are unclear. 

Nevertheless, we believe that in the light of the present day status of the 2(3 decay research 7 
it is useful to discuss and test some "unconventional" (unusual and even strange) ideas and 
approaches to detect the 0z/2/5 decay, like ones presented in this paper, even notwithstanding 
the modest experimental limits reached for the first time and quite uncertain perspectives of 
their advancement from the to-date point of view. 
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